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Abstract 

We demonstrate that in the large N approach developed by the authors in collabo- 
ration with Bardeen, the parameters B & ' and Bg ' parametrizing the K —> 7T7t 
matrix elements (Qe )o and (Q 8)2 of the dominant QCD and electroweak operators 
receive both negative 0(1/N) corrections such that Bq < Bg ' ’ < 1 in agree¬ 
ment with the recent lattice results of the RBC-UKQCD collaboration. We also 
point out that the pattern of the size of the hadronic matrix elements of all QCD 
and electroweak penguin operators Qi contributing to the K —>• ttjt amplitudes Aq 
and A 2 , obtained by this lattice collaboration, provides further support to our large 
N approach. In particular, the lattice result for the matrix element {Qs)o implies 
for the corresponding parameter B g 1 ^ = 1.0 ± 0.2 to be compared with large N 
value = 1-1 A 0.1. We discuss briefly the implications of these findings for 

the ratio e'/e. In fact, with the precise value for Bg ' ; from RBC-UKQCD collab¬ 
oration, our upper bound on B implies e'/e in the SM roughly by a factor of 
two below its experimental value (16.6 ± 2.3) x 10~ 4 . We also briefly comment on 
the parameter Bk and the A I = 1/2 rule. 
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1 Introduction 


The decays K —> irir have played a very important role since their discovery in the 
1950s, both in the construction of the Standard Model (SM) and more recently in 
the tests of its possible extensions. Most of the discussions in the literature centred 
on the following quantities: 


• The ratio 


ReHo 

Re^2 


22.4, 


which expresses the so-called A I = 1/2 rule jTJ[2]. 


( 1 ) 


• The parameter ek, a measure of indirect CP-violation in Kl =► vr7 t decays, 
found to be 

e k = 2.228(11) x 10'V^, (2) 


where cj) £ = 43.51(5)°. 

• The ratio of the direct CP-violation and indirect CP-violation in Kl —> tttt 
decays measured to be [(? ^(ij 


Re(e / /e) = (16.6 ± 2.3) x 10 -4 . 


(3) 


Unfortunately, due to non-perturbative uncertainties originating in the hadronic 
matrix elements of contributing four-quark operators, it took a long time to obtain 
meaningful results for all these observables in QCD. But already in the second half 
of the 1980s, we have developed an approach to K° — K' 0 mixing and non-leptonic 
A'-meson decays [7-lI based on the dual representation of QCD as a theory of 
weakly interacting mesons for large N, where N is the number of colours |12 


15 


The most recent results from our approach can be found in 16 17 
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This approach provided, in particular, first results within QCD for the ampli¬ 
tudes ReAo and Re ^2 in the ballpark of experimental values [lOj. In this manner, 
for the first time, the SM dynamics behind the A I = 1/2 rule has been identi¬ 
fied. In particular, it has been emphasized that at scales 0(1 GeV) long distance 
dynamics in hadronic matrix elements of current-current operators and not QCD- 
penguin operators, as originally proposed in 18 , are dominantly responsible for 


this rule. Moreover, it has been demonstrated analytically why ReAo is enhanced 
and why Re ^2 is suppressed relative to the vacuum insertion approximation (VIA) 
estimates. In this context, we have emphasized that the so-called Fierz terms in the 
latter approach totally misrepresent 1/IV corrections to the strict large N limit for 
these amplitudes 


10 


Our approach, among other applications, allowed us to consistently calculate, 


( 1 / 2 ) 


and 


for the first time within QCD, the non-perturbative parameters Bk, Bq 
B g ' ; governing the corresponding matrix elements of AS = 2 SM current-current 
operator and K — > tttt matrix elements of the dominant QCD-penguin (Qq) and 
electroweak penguin (Qs) operators. These parameters are crucial for the evaluation 
of ek and e'/e within the SM and its various extensions. Other applications of large 
N ideas to K -A 7T7r and Bk, but in a different spirit than our original approach, 
are reviewed in 119]. We will comment in Section 6 on those which reached very 
different conclusions from ours. 

It is interesting and encouraging that most of our results have been confirmed by 
several recent lattice QCD calculations that we will specify below. While the lattice 
QCD approach has a better control over the errors than our approach, it does 
not provide the physical picture of the dynamics behind the obtained numerical 
results. This is in particular seen in the case of the A I = 1/2 rule where our 
analytic approach offers a very simple picture of the dynamics behind this rule, as 
summarized again in [16,17 . 


In the present paper, we briefly compare in Section 2] the status of lattice results 
for Bk and the A/ = 1/2 rule with the ones obtained in our approach. Subse¬ 
quently, in Section 3] we demonstrate that the pattern of the size of the matrix 
elements for penguin operators presented recently by the RBC-UKQCD collabora- 
and A 2 amplitudes 21 gives another support to our approach. In 


tion for Aq 20 


0 (1/2) (3/2) 

we derive upper bounds on the parameters Bq ' ' and B g ' ' and discuss 

briefly in Section 5j their phenomenological implications for s' /e. In Section 6: we 

describe briefly the results obtained in other large N QCD approaches. An outlook 

is presented in Section 7| 


2 Bk and the A I = 1/2 Rule 
2.1 13 k 

The scale and renormalization scheme dependent parameter Bk(^) is related to the 
relevant hadronic matrix element of the AS = 2 operator 


Q = (sd)v-A(sd)v-A 


(4) 
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as follows!) 

(K°\Q^)\K°) = B K {^-F 2 K m 2 K . (5) 


More useful is the renormalization group invariant parameter Bk that is given 

by (22) 


Bk = B k {h) 



—d 

«i 3 ) (h) 



4lT 


J3 , 


9 (IV — 1) 
N(11N — 6 ) ' 


( 6 ) 


We have shown the ^-dependence of the exponent d in the leading term to signal 
that d vanishes in the large N limit. The coefficient J 3 is renormalization scheme 
dependent. This dependence cancels the one of Bk( fJ-)- 

As in the strict large N limit the exponent in ( 6 ) and the NLO term involving 
J 3 vanish, one finds (7] that independently of any renormalization scale or renor¬ 
malization scheme for the operator Q 


Bk —> 0.75, (in large N limit, 1986). 

It can be shown that including 1/N corrections suppresses Bk so that 
Bk < 0.75 , (in 1/N expansion). 

Our latest analysis in our approach gave [16] 


23 


(7) 

( 8 ) 


B k = 0.73 ± 0 . 02 , (in dual QCD), 


(9) 


where the error should not be considered as a standard deviation. Rather, this 
result represents the range for Bk we expect in our approach after the inclusion of 
NLO QCD corrections and the contributions of pseudoscalar and vector mesons as 
discussed in detail in 


16 


On the other hand, the world lattice average for Bk based on the calculations 

1 calculations (recent FLAG update 


of various groups 24 


of 27 


29 


reads for Nf = 2 
B k = 0.766 ± 0.010, (in lattice QCD, 2014). 


( 10 ) 


See also the recent analyses in |;30,||31 . While this result violates the bound in ( 8 j), 
it should be noted that a number of lattice groups among [24 29 published results 
with central values satisfying the bound in ( 8 ) but the errors did not allow for a 
clear cut conclusion. In fact, the most recent update from staggered quarks |31 


quotes precisely Bk = 0.738 ± 0.005 but additional systematic error of 0.037 does 
not allow for definite conclusions. Similarly, the Rome group 32 finds basically 
the result in (9). We expect therefore that improved lattice calculations will satisfy 
our bound one day and in a few years from now lattice average for Bk will read 
B k ~ 0.74. 


Finally, let us remark that while the lattice approach did not provide the expla¬ 
nation why Bk is so close to its large N limit 0.75, in our approach the smallness of 
1/N corrections follows from the approximate cancellation of negative pseudoscalar 
meson contributions by the positive vector meson contributions. 


1 In this paper, we use the normalization of weak decay constants given in (31). 
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2.2 A/= 1/2 Rule 


A very detailed comparison of the calculations of Redo and Red 2 in our approach 


and the lattice QCD has been presented in 16 and our present discussion is meant 


to be an update due to new results of the RBC-UKQCD collaboration on Redo |20j . 

First, let us mention that both the dual approach to QCD and lattice approach 
obtain satisfactory results for the amplitude Red 2 . On the other hand, whereas we 
find 116 

/ Redp \ 

VRedJ, 


= 16.0 ± 1.5, 

dual QCD 

the most recent result from the RBC-UKQCD collaboration reads |20 


(ID 


Redo\ 
Red 2 ) 


= 31.0 ±6.6. 


( 12 ) 


lattice QCD 


Due to large error in the lattice result, both results are compatible with each 
other and both signal that this rule follows dominantly from the QCD dynamics 
related to current-current operators. But our approach, being analytic, allows to 
connect the A I =1/2 rule to the main properties of QCD: asymptotic freedom and 
the related evolutions of weak matrix elements which at long distance scales can 
be performed in the dual representation of QCD as a theory of weakly interacting 
mesons for large N. As lattice QCD calculations are performed basically at a 
single energy scale, no such physical explanation of this rule is expected from that 
framework. To this end, lattice calculations would have to be performed at scales 
below 1 GeV which is straightforward in our approach but appears impossible by 
lattice methods at present. 

On the other hand, from the present perspective only lattice simulations can 
provide precise value of Redo one day, so that we will know whether some part of 
this rule at the level of (20 — 30)%, as signalled by the result in ([TTJ) , originates in 
new physics (NP) contributions. Indeed, as demonstrated in [33] , a heavy Z' and in 
particular a heavy G' in the reach of the LHC could be responsible for the missing 
piece in Redo in (11}. On the basis of the analysis in [33] it is much harder to 
bring this ratio with the help of NP from 31 down to 22 without violating A Mk 
constraint, but this requires a separate study. 

Of some interest is the ratio of the matrix elements (<^ 2)0 and (Qi)o- It equals 
—2 in the large N limit, corresponding to /i = 0 16 . Evolving these matrix elements 


to /r = 1 GeV in the meson theory and subsequently to /1 = 1.53 GeV in the quark 
theory, we find in the NDR-MS scheme 2 

(Q2)o 


(Q 


= -1.50 ±0.10, 


1/0 


H = 1.53GeV, (dual QCD). 


The corresponding result in [20] reads 

(Q2)o 


(Q 


1/0 


= —1.12 ±0.49, n = 1.53 GeV, (lattice QCD). 


(13) 


(14) 


In view of large uncertainty in the lattice result, these two ratios are compatible 
with each other. We expect on the basis of the results in (llj) and (12) that this 
ratio will be eventually found in the ballpark of —1.4. 


2 We thank Martin Gorbahn for checking this result. 
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3 Matrix Elements of Penguin Operators 

3.1 Preliminaries 

We will consider the usual basis of operators contributing to K —>• 7 T 7 t amplitudes 


34 , namely 

Current—Current: 


Q i = ( s a up)v~A ( upd a )v-A , Q2 = ( su)v-a ( ud)v-A 


QCD—Penguins: 

Q 3 = ( sd)y-A 22 (qq)v-A, 

q=u,d,s 


Qa = ( S a dp)v-A '22 ( Q/3Qa)v-A 

q=u,d,s 


(15) 


(16) 


Qb = ( sd)v-A 22 (mW+a, Qe = ( s a dp)v-A 22 (QPQ<*)v+a (17) 


q=u,d,s 

Electroweak Penguins: 


q=u,d,s 


O 7 = X ( sd)v-A 22 e 4 (w)v+A, Qs = X ( s a dp)v-A 22 e q(W<*)v+A (18) 


q=u,d,s 

Q 9 = ^ ( sd)v-A 22 

q=u,d,s 


q=u,d,s 


e q {qq)v-A, Q 10 = ^ ( S adp)v-A 22 

q=u,d,s 


eq (qpqa)v-A (19) 


Here, a, f3 denote colour indices and e q denotes the electric quark charges reflecting 
the electroweak origin of Q 7 ,..., Qio- Finally, (sd)v-A = s«7m(1 ~l5)d a as in 
Recently, the RBC-UKQCD collaboration published their results for the matrix 


elements (Qi)o 20 . Their matrix elements are given for three dynamical quarks at 


p = 1.53 GeV, which is too high for the direct comparison with our approach in 

the case of current-current operators. On the other hand, the parameters B g 1(/ ^ 
( 3 / 2 ) 

and FI ; of the QCD penguin operator Q& and the electroweak penguin operator 


Qs are known 34 to be practically scale independent for 1.0 GeV < p < 3.0 GeV. 


Therefore these results constitute a useful test of our approach. Another issue is 
the colour suppression of some matrix elements of other penguin operators which is 
predicted within our approach. We would like to check whether the pattern of this 
suppression is also seen in the lattice data. 


3.2 Hadronic matrix elements 

The hadronic matrix elements of operators Qi that are most useful for our discus¬ 
sions are 

(Qi)i = ((tttt)/ \Qi\ K) , (20) 

with 7 = 0,2 being strong isospin. 

It should be recalled that for /x < m c , when charm quark has been integrated out, 
only seven of the operators listed above are independent of each other. Eliminating 
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then Q 4 , Qg and Q 10 in terms of the remaining seven operators allows, in the isospin 
symmetry limit, to find the following important relations 34p 


(Q4)o — (Q3)o + (< 22)0 — (<2i)o ) 
((29)o 

(Qio)o 

(09)2 


^(^1)0 “ 2 ^ 3 )° > 

(<22)0 + ^{ Qi)o - 2(* 2 3 )° > 


where we have used 


(Qio )2 — 


{ Qi )2 = (<22)2 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


We have checked that these relations have been used in [20]. 

Of particular importance for our discussion are the matrix elements 


(< 2 6 ( m))o h 

{Q&(H ))2 = 


2 




2 m|- 


(F k - F n ) b£ /2) 

■ ( 3 / 2 ) 


2\/2 [m s (n) + 

2 


F 

r 7T -Dg 


(Q 8(p))o — 7T 


2ro| 


2 [m s (p) + m d (n)_ 


p pi 1 / 2 ) 
r 7 r -°8 


with 7 8,35 


Bj 1/2) = ^ 3/2) = S^ 1/2) = 1, (large N Limit). 


(26) 

(27) 

(28) 

(29) 


Note, that using the definition of Bi parameters consistent with the large N limit 
of QCD, as given above, implies that their values in the VIA |34| read 


B , 


( 1 / 2 ) _ 


Bf /2) sa 0.99, 


6 — 1, -^8 ~ ^8 


B, 


( 1 / 2 ) 


1.2 


(VIA). 


(30) 


We will return to this point in the next section. 


The input values of parameters entering these expressions are given by 136,37 


Fk 

iV = 130.41(20) MeV, -//- = 1.194(5) 

FK 

m s (m c ) = 109.1(2.8) MeV, m d (m c ) = 5.44(19) MeV . 


(31) 

(32) 


It should be emphasized that the overall factor h in these expressions depends 
on the normalisation of the amplitudes Aog. In [34] and recent papers of the RBC- 
UKQCD collaboration j; 2lj|38| h = ^/3/2 is used whereas in most recent phenomeno¬ 
logical papers 16,19 33 39 , h = 1. In the present paper we will keep general h so 
that, e.g., the decay amplitude K + reads ( 3 / 2 h)A 2 - 


3 In writing (21) we neglect a small 0(a s ) correction in the NDR scheme which is explicitly given in 
(4.44) of “ 


34 
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Comparing the expressions (26) and (28) with the lattice results in [20], we find 
(see also 1 40 ) 4 


B, 


( 1 / 2 ) 


= 0.57 ±0.19, 


B, 


( 1 / 2 ) 


= 1.0 ± 0.2 , (lattice QCD). (33) 


On the other hand, comparing (27) with the value for this matrix element obtained 


(34) 


by RBC-UKQCD collaboration in [21 one extracts [39] 

B^ 2 ' 1 = 0.76 ± 0.05 , (lattice QCD). 


All these results are very weakly dependent on the renormalization scale. The 
quoted values correspond to p = 1.53 GeV. Basically, identical results are obtained 
for p = m c used in [40[. However, as stated before (33), in extracting these param¬ 


eters from [ 20 ] it is important to use the quark masses at that scale. 

As we will demonstrate in the next section, these lattice results are consistent 
with the large N approach. Indeed, we will show that the following pattern emerges 
at next-to-leading order in our dual approach: 

Fir 


b£ /2) = 1 - 


Fk ~ F n 


O(-) 

k n j 


< 1 , 


B 


( 3 / 2 ) 


B 


( 1 / 2 ) 


= l- 0 (-)cl. 

= 1 +0(i) >1 


(35) 

(36) 

(37) 


( 3 / 2 ) 

We would like to recall that strong indication for the suppression of B$ ' below 

while in the case 


unity in our approach have been found already in 1998 in 41 

of B g 1//2 ^ no clear cut conclusions could be reached. Our present analysis of both 
( 1 / 2 ) ( 3 / 2 ) 

Bq 1 ’ and Bg 7 7 clearly indicates the negative signs of 1/N corrections to the 
leading result in (29). 


Finally, the lattice results in 20 and 1211 exhibit colour suppression of the 


matrix elements of Q 3 , Q 5 and Qj operators relative to the ones of Q 4 , Q§ and Qs, 
respectively: 

'Fk 
F K 

{Ql)2 


(Qs)c 

(Q4)c 


= -0.18 ±0.25, -^-1 


(Q 7)c 
(Q 8)c 


= 0.13 ±0.04 


(Qs)q 


(05)0 

^5/0 = 0.10 ±0.05, 
(<36)o 

(38) 

= 0.22 ±0.01. 

(39) 


These results are consistent with the large N approach. Indeed, as we will demon¬ 
strate soon, the ratios in (38) are 0(1/N 2 ) while the ratios in (39) are 0(\/N). 

These results allow to simplify some of the relations between the matrix elements 
so that it is justified to use the relations 

(Q 4 )o = (Q 2 )o — (<3l)o ; (40) 


(139)0 — ^( Q 1 )o > 

(< 3 lo)o = (132)0 ± ^(Ql)*)’ 

which simplify the phenomenological analysis of e 7 /e in [401. 


(41) 

(42) 


4 To this end, the values m s = 102.27 MeV and rrid = 5.10 MeV at /r = 1.53 GeV have to be used. 
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4 Derivations 


The large N numerical values of the |A5| = 1 matrix elements already displayed 


in Section 3.2 are most easily derived from the effective theory for the the pseudo- 
Goldstone field 

U(tt) = exp (*-\/2y) (43) 

with 7r = A a 7r a , the meson nonet lying below the one GeV and /, the associated weak 
decay constant scaling like y/~N. In particular, the electroweak penguin operator 
introduced in (18) and “Fierzed” into a product of two colour-singlet quark densities, 
namely, 

Qs = -12 V' ( SLqR)e q (qRd,L ) (44) 


q=u,d,s 


can be hadronized by considering the leading chiral effective Lagrangian in the large 
N limit: 


f 2 r 

L e f{ (p 2 , N) = J —Tr d li Ud fl U + + r(mf/ t + Urn)) 


(45) 


Indeed, a straightforward identification of the second term in this equation with the 
standard Dirac mass term in QCD 


TqcD(mass) = ~{qLmq R + qRm ] q L ) 

(46) 

allows us to hadronize all colour-singlet quark densities 



(47) 

?LlR = -LrU' b “ 

(48) 

such that 

Qs = -^/V U d «e q U^ s . 

q=u,d,s 

Consequently, the factorized matrix elements of the A I = 
components of Q$ in the large N limit are 

(49) 

1/2 and A/ = 3/2 

(Qs)o = \fr 2 

(50) 


(51) 


Similarly, the QCD penguin operator Q§ introduced in (17} and “Fierzed” into 
a product of two colour-singlet densities reads 

Q 6 = -8 ( s L q R )(q R d L ) = -^/V U dq U^ qs = 0 . (52) 

q=u,d,s q 

As a matter of fact, one has the relation 


r(n) = 


2 m 2 K 


m s (n ) + md(n) 


(53) 
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at the level of L e s(p 2 ,N). Yet, at this level, the absence of SU( 3) splitting among 


the weak decay constants implies ill-defined (Q s)o ,2 matrix elements in (50) and 


(51) as well as a vanishing Qq operator in (52). 


It is well known [8,42 that the next-to-leading term in the chiral effective La 
grangian 

4 f 2 r 




md 2 U t + d 2 Um) 

solves both problems since it leads to realistic weak decay constants 


(54) 


Ft 


= (! + ^)/ 
x 


F k 

/ = 1 + 

-L 7T 


2 2 
nix - m% 

V 


thanks to its derivative dependence and, simultaneously, it implies 

Q °=-T{i) 2(Bf ‘ ua * u ' )d ‘ +0( -k ) 

from the shift induced in the hadronized quark densities 


-1 yj 


U 


u -^d 2 U 

A l 


(55) 


(56) 


(57) 


(58) 


through its mass dependence. Taking these corrections into account, we now repro¬ 
duce the large N matrix elements given in 


(26) 


with the normalization (29) 


for the i?6,8 coefficients if contributions 0(?n 2 /A 2 ) to Qg are neglected 

At this point, it is worth emphasizing that, here, we consistently normalize the 
|AS| = 1 Bqs to unity in the large N limit. Such is unfortunately not the case 
for the |AS| = 2 Bk parameter conventionally normalized to one with respect to 
VIA in Q. Had the AS = 2 matrix element in (5) been normalized relative to its 
large N value, the most precise Bk parameter extracted from lattice QCD in (10) 
would read Bk = 1-021 ± 0.013 nowadays and our result in (9) Bk = 0.97 ± 0.03. 
In 


34,41 , Bq and Bg were also normalized with respect to the VIA as in (30). 


We are now in an ideal position to estimate 1/N corrections encoded in the 
Bq^ parameters. The factorizable 1/N corrections to |A S\ = 1 density-density 


operators are fully included [35] in the running of quark masses in (53). Let us thus 
focus on non-factorizable one loop corrections induced by L e g(p 2 , N ). Applying the 


background field method of 43 , we find 

16 ln(A 2 /M 2 ) 


U d m Ws (A) = U dq U Ws (M ) - 


/ 4 (4vr/) 2 


2J ds ji q + {J L J L ) ds d q,q ] (M) (59) 


with A = 0(1 GeV) the euclidean ultraviolet cut-off of the effective theory (45) to 


be matched with the non-factorizable short distance evolution, M = 0(rriK ) and 


*<& = it(d„UU'r b 

(60) 

fa A = <F(afaur b 

(61) 
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the colour-singlet left-handed and right-handed hadronic currents derived from 
Lettip 2 j IV), respectively. 

Applied to the specific K —> -rnr decay processes, 


the first (L-R) current-current operator in (59), also present with the right 


relative sign in the VIA through a Fierz transformation, does not contribute 
to the matrix element (Qe)o since 


V3 

2 


Tr(J R ) = 


(62) 


the second (L-L) current-current operator in (59), absent in the VIA, does not 
contribute to the matrix elements (Qs)o ,2 since 


Tr(e g ) = 0. 


(63) 


An explicit calculation of the surviving Q 4 and Q 7 matrix elements (in the large 
N limit) gives then, respectively, 


b< 1/2 > = 1 - - 
6 2 


4 1/2) = 1 + 


B {3/2) = 1-2 


Fk - F T[ 

(m? K - ml) 
(4vrE, 

( m K ~ ml) 


mi mi 


A 2 A 2 

2 ln(l -|—^ 2 ) = 1 + 0.08 ln(l + 


m\ 


8 


m\ 


8 


(47TF; 


A 2 A 2 

ln(l + -) = !- 0.17 ln(l + 


m s 


ms 


(65) 

( 66 ) 


with pseudoscalar mass scale parameters bounded necessarily by the effective cut-off 
around 1 GeV: 

me, 8<A. (67) 

First, we emphasize most important properties of these results: 

• For A = 0, corresponding to strict large N limit and matrix elements evaluated 
at zero momentum, B^"' = B^ 2 ' 1 = B = 1 in accordance with (29). 

( 1 / 2 ) ( 3 / 2 ) 

• With increasing A, the parameters Bq ' ; and B g ' decrease below unity and 

( 1 / 2 ) ( 3 / 2 ) 

Bq 7 7 decreases faster than B$ . Consequently, at scales 0(1 GeV) relevant 

( 1 / 2 ) ( 3 / 2 ) 

for the phenomenology both B g ' and B g ' 7 are predicted to be below unity 

( 1 / 2 ) ( 3 / 2 ) 

and there is strong indication that B§ < B% . 

( 1 / 2 ) ( 3 / 2 ) 

• While the dependence of Bq ' 7 and B g “ on A < 1 GeV is stronger than their 
dependence on /r in the perturbative regime, these two properties of B^ 2 ' 1 
and .Bg ' 7 are at the qualitative level consistent with the numerical analysis 
performed for B g 1 ^ and B^ 2 ^ by means of the standard renormalization 


group running m 


34 


Indeed as seen in Figs. 11 and 12 of that paper B { 


( 1 / 2 ) 


( 3 / 2 ) 

decreases with increasing /r, faster than B & ' 7 , albeit in this perturbative range 
the dependence of Bg 7 7 and Bg ' 7 on fi is very weak. While the analysis 
in (341 includes NLO QCD and QED corrections, the inspection of the one- 
loop anomalous dimension matrix allows to see these properties explicitly. In 
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particular Qq mixes with the linear combination (Q 4 + QQ and we find for 
Ail < fJ-2 < m c 


4 1 / 2 , te) = 4 1/2, («) 


2 vr Vi 


1 + 


(Q4(Aij))o 
{Q&(^l))o ) _ 


( 68 ) 


From (57) |(<^6 (a^i)) oI > |(<^ 4 (/^l))oI such that B g 1 ^ decreases with increasing 
/x. On the other hand, in the LO Qs runs only by itself and the one-loop 
anomalous dimension matrix implies 


£i I/2 , 3/2, ta) = B? /2 ' 3/2 Vi), 


(69) 


which follows from exact SU(3) symmetry imposed in SD calculations. The 

breakdown of SU(3) is only felt in the matrix elements of Q § making in the 
( 3 / 2 ) 

LD range B% ~ dependent weakly on the scales involved. In view of this, the 
( 1 / 2 ) ( 3 / 2 ) 

suppression of both Bg 7 and B% 7 7 below the unity can be considered as a 

solid result and our explicit calculation as well as different behaviour of Qq 

( 1 / 2 ) ( 3 / 2 ) 

and Qs under flavour SU(3) provide a strong support for B^ 7 < B$ . On 

the other hand, 


B, 


( 1 / 2 ) 


, 8 ~ [^ 3/2) ]~ 1 /2 ( 70 ) 

slightly increases with A which is also consistent with the standard renormal¬ 
ization group running [40] , 

Next, we observe that: 

• The numerical value of the parameter B g 1 ^ suffers from rather large uncer¬ 
tainties. This feature is related to the fact that Qq vanishes at leading order in 
chiral perturbation theory (see (52)). The 1/IV logarithmic correction in (64) 
is therefore artificially enhanced by the factor F n /(FK — F n ) « 5 such that 


B, 


( 1 / 2 ) 


< 0 . 6 . 


(71) 


The parameter B^^ has a very small 1/IV correction. At 0{\/N 2 ), one larger 
contribution might arise from the anomalous effective Lagrangian 


f 


L eS (p °, 1/IV) = J — ( ^ ) Tr(ln U — lnt/t) 


32 V N 


i 2 


(72) 


that solves the so-called U(1)a problem 44 by providing the 7 / pseudoscalar 

(73) 


with a physical mass in the large N limit 45 


2 1 2 o 2 

m / +m„ — lm K 


m,n 


0.7 GeV 2 


Applying again the background field method, we obtain 


U dq U Ws (A) = 


jV(4vr/) 2 y AP } 


Such a negative contribution to Qs has been included in 41 
consistent estimate beyond 


B, 


( 1 / 2 ) 


1 


U dq U Ws (M). (74) 

However, any 
(75) 


would require a full calculation at 0(1/N 2 ). 
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The parameter B^ 2 \ for which the l/N expansion is more reliable, is found 
in the range 

(76) 


0.7 < B, 


(3/2) 


-^8 — 


< 0.9 


if 7718 > 777 K- 


Lattice result for B^ 2 ^ in (33) turns out to almost saturate our bound. But one 


should realize that although we are confident about the suppression of B g 1 ^ below 
unity, its actual size is rather uncertain. For instance the inclusion of dynamical 
ars presently frozen in Ay could reduce the coefficient in front of the logarithm 


sea. 


in (64) making B ^ 2 ' 1 larger. This uncertainty in the value of B g 1 ^ explains also 
why it took so long to calculate B g 1 ^ in lattice QCD even with a large uncertainty 

as seen in (33). On the other hand, the range for B% in (76) is consistent with the 

(1/2) (3/2) 

(34). These results indicate that indeed B1 7 7 could be smaller than B1 7 . 


one m 


Yet, in view of the large numerical uncertainties in the case of B^ 2 \ we cannot 
(1/2) (3/2) 

exclude that Bq 7 7 is as large as B g 7 . We therefore believe that the best way 
of summarizing our results for B g 1 ^ and B g 3 ^ is given in (35) and (36) together 
with 

" < Bf /2) < 1. (77) 


B, 


( 1 / 2 ) 


6 — 


Below 1 GeV we have seen in (59) that density-density operators transmute into 


current-current ones at 0(1/N). But power counting in our effective theory does 
not allow the other way around, namely current-current operators evolving into 
density-density ones. This is fully consistent with the evolution of hadronic matrix 


elements above /j, = 1 GeV studied already in 34 and is opposite to the evolution of 


the corresponding Wilson coefficients. Now, in the large N approach, it has already 
been shown flu, 43j| that any (L-L) current-current operator evolves as 


1 


Jf Jf(A) = Jf Jf (0) - O(-) [ 2 Jf Jf - 5 ad (J L J L ) cb - 5 cb (J L J L ) 
to stand in contrast with the wrong relative sign in the VIA analogue 


ad 


l 


Jf Jf (A) = Jf Jf (0) + - Jf Jf (0) 


(0) (78) 


(79) 


As a consequence, summing over c = d = u, d, s we conclude that the matrix element 
{Qs)o which vanishes in the large N limit is formally 0(l/N 2 ) relative to (QQo in 
our effective theory. Such a strong suppression could have been anticipated from 
the LO short-distance evolution of the four-quark operator Q 3 into another linear 
combination of Q 4 and Qq: 


Qz{H 2 ) — <33(^1) - Cf —) 


2 Q4 + Qa 


(Mi)- 


(80) 


At long distance, the further Q 3 evolution undergoes an important numerical cancel¬ 
lation since (57) tells us that the Q 4 and Qq operators are not independent anymore. 
Following (|78), this numerical cancellation is not a mere coincidence but the result 
of a consistent l/N expansion. 


In the same manner, it has been proved 43 that any (L-R) current-current 
operator evolves as 


«(*) = Jt b Jg(0) + 0 ( 4 ) 


U ad { 8 tf) cb + 8 U ad {U^) cb 


( 0 ) 


(81) 
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with 5U proportional to (□[/ — UOU^U). 

Consequently, the matrix element (Q 5)0 which vanishes in the large N limit is 
0(1/N 2 ) relative to F n /(Fx — F- k )(Qq)o in our effective theory. As already men¬ 
tioned, an enhancement factor has to be introduced to compensate for the “acci¬ 


dental” chiral suppression of (Qq )0 in (52). On the other hand, the matrix elements 


(C?7>o = ^Fv{m 2 K -ml) 

{Qt )2 = ~-^=F n (m 2 K - ml) 


(82) 


(83) 


are 0{p 2 ) but non zero in the large N limit. With the matrix elements (<3s)o,2 given 
in (50), (51) and being O(p 0 ), the resulting ratios (< 37 ) 0 , 2 /(<3s)o,2 are at the level 


of a few percent and can thus be neglected. The long-distance evolution of Qj in 
(81) leads then to matrix elements proportional to {Qs)o, 2 , though 0(p 2 ). This is 


clearly at variance with its LO short-distance evolution, namely 

< 37 (^ 2 ) = Q7(^i) + <9( —) <38 (at) • 


(84) 


As already explicitly stated in [43], this suggests the necessity to introduce higher 
resonances beyond our effective theory truncated to the low-lying pseudoscalars. In 
a dual representation of QCD the matrix elements (< 37 ) 0,2 should then be dominantly 
O(p 0 ), but 1/AA-suppressed, with the bound 


(Q 


7/0 


(Q 


< V2 


7/2 


(85) 


resulting from the isospin decompositions in (82)-(83) and (50)-(51). 


5 Implications for s'/s 


We will now briefly discuss the implications of our results for e'/e. To this end 
we will use the analytic formula for s '/e in the SM derived recently in |40]. In 
obtaining this formula it has been assumed that the SM describes exactly the data 
on CP-conserving K —> nir amplitudes: ReAo and ReA 2 . This allowed to determine 
the contributions of the (V — A) tg> (V — A) QCD penguin operator <3 4 and of the 
electroweak penguin operators <39 and < 3 io to e'/e much more precisely than it is 
presently possible by lattice QCD and large N approach. This determination was 
facilitated by our results on the suppression of the matrix element (< 33)0 implying 
the relations (40)- (42) 5 . 

The formula in question reads [40] 


- = 1(T 4 

£ 


ImAt 


1.4- 10 


-4 


o(l - n eS ) ( - 4.1 + 24.7 Rg 1/2) ) + 1.2 - 10.4 B, 


>(3/2) 


where 


ImAt = Im(' V td V* s ) = \V ub \\V cb \ sin 7 


( 86 ) 

(87) 


The final numerical analysis in 40 leading to ( 86 ) included also small corrections from Q 3 and other 


corrections from subleading operators. 
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and 140 46-48 


a = 1.017, H e ff = (14.8 ± 8.0) X 1(T 2 . ( 88 ) 

\V u b\ and | V c b | are the elements of the CKM matrix and 7 is an angle in the unitarity 
triangle. The parameters a and represent isospin breaking corrections 
See these papers and [40 for details. 


46 


48 


(1/2) (3/2) . 

Setting all parameters, except for IT ' and Bo , in ( 86 ) to their central 


values we find 


Re(e / /e) = 8.6 x 10“ 4 , 

Re(e / / e ) = 6.4 x 10 -4 , 
Re(e , /e) = 2.2 x 1 (T 4 , 


(B^ /2) = 1.0 ,Bf /2) = 1.0) 
(b£ /2) = 0.8, Bf /2) = 0.8) 
(B^ /2) = 0.6, Bf /2) = 0.8). 


(89) 

(90) 

(91) 


A detailed anatomy of s'/e in the SM is presented in [40j, where various uncertainties 

related to NNLO QCD corrections and other uncertainties are discussed. But these 

three examples indicate that taking our bounds into account and guided by the 
(1/2) (3/2) 

results on and B\ ' from lattice QCD and our dual approach, our SM 

prediction for s'/s appears to be significantly below the data given in <©• 

In each of these predictions, there are uncertainties from the value of ImAt, the 
unknown complete NNLO corrections to Wilson coefficients of contributing opera¬ 
tors, a s , mt and other input parameters. But they appear not to change the conclu¬ 
sion that, presently, the SM prediction for s'/s is significantly below the data. Our 
upper bound on B g 1 ^ plays an important role in this result as otherwise increasing 
B g 1 ^ above unity would allow to fit easily the data. 


6 Comments on other large N QCD approaches 

In [49] the authors analyse A I = 1/2 rule and s'/s in the chiral limit including 1/N 
corrections. Their results differ drastically from our results. In particular, their low 
energy ’’Extended Nambu-Jona-Lasinio” model (ENJL) gives 

b£ /2) » 3, B ( 8 3/2) « 1.3, (ENJL) (92) 

namely B^ 2 ^ roughly by a factor of five larger than lattice calculations and our 
results. With such high values of B^ 2 \ QCD penguins play an important role in 
the explanation of the A I = 1/2 rule and the experimental value of e'/e can easily 
be reproduced, again in contrast with lattice QCD and our dual QCD approach. 

In |50| the authors consider a low energy model including the light pseudo-scalar, 
vector and scalar poles only in the chiral limit. Within this so-called ’’Minimal 
Hadronic Approximation” (MHA), they also obtain much larger values of B g 1 ^ 
and B^ 2 ' 1 than in our approach 

B^ ,2) fn 3, B {3/2) fn 3.5, (MHA) (93) 

and find then good agreement with data for s'/s. 

In [51fj53j the authors rely on dispersion relations and ’’Finite Energy Sum 
Rules” (FESR) in the chiral limit to extract the electroweak penguin matrix el¬ 
ements from ALEPH and OPAL data. Doing so, they obtain central values for 
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(3/2) 

Rg ' ' by a factor of two to three larger than in lattice QCD and our dual QCD 
approach 


1.3 < i?g 3/2) < 2.5, (FESR). 


(94) 


We conclude that the recent lattice QCD results tend to demonstrate that the 
various models considered in 49 - 53 do not represent properly the low energy QCD 
dynamics at work for penguin matrix elements, but confirm the structure of our 
dual QCD approach. 


7 Summary and outlook 


In the present paper, we have compared the structure of the hadronic matrix ele¬ 
ments in K —» 7T7T decays obtained within the dual approach to QCD with the one 
obtained recently by the RBC-UKQCD lattice approach to QCD and commented 
briefly on the status of the parameter Bk and the A I = 1/2 rule. Our main results 
are as follows: 


• The status of Bk is very good as both our approach and the lattice QCD 
calculations give this parameter very close to 0.75. But we expect that Bk 
from the lattice approach will decrease by a few % in the coming years. 


• While the results for Re %2 obtained in both approaches agree well with the 
data, the central value of ReAo from RBC-UKQCD collaboration is by a factor 
ot two larger than in our approach and 40% above the data. While our result 
in © appears from present perspective to be final in our approach, significant 
improvement on the lattice result is expected in the coming years. This will 
allow to find out whether at some level of 20 % new physics could still be 
responsible for the A I = 1/2 rule. An analysis anticipating such possibility 
has been presented in [33]. 


(3/2) 

As the upper bound on B g 7 in (36) has been already indicated in 


41 


(1/27 


one 

Our 


of the most important results of our paper is the upper bound on B g 

(1/2) (3/2) 

estimate suggests that Bq ’ < Rg ' < 1 , but the precise values can only be 

obtained by lattice methods. 


• Among other results of our approach supported by recent results from RBC- 
UKQCD is the strong suppression of (<33,5(/-0) and Rg 1 ^ 2 ^ ~ 1 . 

If indeed the emerging pattern B g 1 ^ < Rg 3/ ^ < 1 with Rg 3// ^ = 0.8 ±0.1 will 
be confirmed by more precise calculations one day, the very recent analysis in 1 40 1 
and our paper show that e'je within the SM will be found roughly by a factor 
of two below the data. For a detailed phenomenological discussion of the state 
of e'Ie within the SM including all errors and future theoretical and experimental 
prospects we refer to |40j. On the other hand, first phenomenological implications 


of our results on new physics models have been presented in 54,55 
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